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We perform a simplified analysis of finite temperature properties of the generalized Falicov- 
Kimball model witli Hund coupling. The model is studied on the square lattice for fixed values 
of f-electron and d-electron densities p/ = 1/2, pd = 1. Using an approxiamte scheme we estimated 
free energies of the following three phases: the charge and spin ordered (SCO), the charge ordered 
(CO) and the non-ordered (NO). Comparing the free energies we detected phase transitions and 
found how the transition temperatures depend on the on-site interaction parameters, ft appears, 
that the transition temperture between SCO and CO is much lower than between CO and NO, but 
their maximum values are attained at values of coupling constants very close to each other. 

PACS numbers: 71.10.Fd, 71.28.4-d, 73.21.Cd, 75.10.-b 



I. INTRODUCTION 



Charge nad magnetic order are observed in many correlated electron systems as, for example, in R^^xSr^NiOA, 
where R = La,Nd One of the models capable to describe such superstructures is the Falicov-Kimball model, 
extended by the spin-dependent local interaction, that reflects the first Hund's rule 0. 

Until now the model has been studied only at zero temperature. Using the method of restricted phase diagrams 
or small cluster numerical diagonalization it was shown how the charge distribution and spin arrangement change 
with parameters of the local interaction and densities of localized and itinerant electrons 0, Q. In particular, for 
Pf = 1/2 and = 1 it was found that in 2D the ground state forms the checkerboard charge pattern consistent with 
the simplest antiferromagnetic order [3| . 

Here we investigate behaviour of the system at finite temperatures. Our purpose is to determine the way of 
transformation of the ordered phase into disordered one when temperature increases. The transformation appears to 
split into two phase transitions: at Tmo between the low temperature phase, where both charge and spin are ordered 
(SCO) and the intermediate phase, with only charge but not spin ordered (CO), and then at Tco between CO and 
the high temperature phase, that is disordered with respect to both the spin and charge (NO). A visualization of 
arrangements of the localized electrons in the three phases is shown in Fig. 1. 
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FIG. 1: Typical arrangements of the localized electrons (snapshots) in three different phases: a) the spin and charge ordered 
(SCO), b) the charge-ordered (CO) and c) the non-ordered (NO). 

The model Hamiltonian is 

-JE«T-<i)KT-"a)' W 

i 

where < i,j > means the nearest neighbor lattice sites i and j, a and rj are spin indices, di^cr {d^a) annihilation 
(creation) operator, and nf^ ('^f (j) is an occupation number of itinerant (localized) electrons. The on-site interaction 
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between localized and itinerant electrons is represented by two coupling constants: (7, which is spin-independent 
Coulomb-type and J, which is spin-dependent and reflects the Hund's rule force. The hopping amplitude t we set 
equal to one, so we measure all energies in units of t. 

Double occupancy of the localized electrons is forbidden, implying the on-site Coulomb repulsion Uf / between two 
f-electrons is infinite. Consequently, at a given site the f-electron occupancy is assumed to he Uf = 't-/,! + '^z,! ^ 1 
and the d — electron occupancy to be Ud — nd.-\ + nd,i < 2. So there are 3 states per site allowed for the f-electrons 
(uf = 0; riy i = 1 and rif^i = 0; nj^| — and ny.j = 1) and 4 states per site allowed for the d-electrons (rid — 0; 
Ud.^ = 1 and Ud.x = 0; rid,]- = and Ud.i = 1; = 2). 

All single-ion interactions included in U) preserve states of the localized electrons, i.e. the itinerant electrons 
traveling through the lattice change neither occupation numbers nor spins of the localized ones. In other words, 
[H, f^firj] ~ for all i and 77, so the local occupation number is unchanged. 

The localized electrons play the role of an external, charge and spin dependent potential for the itinerant electrons. 
This external potential is "adjusted" by annealing, so the total energy of the system attains its minimum. In other 
words, there is a feedback between the subsystems of localized and itinerant electrons, and this is the feedback that 
is responsible for the long-range ordered arrangements of the localized ones, and consequently for the formation of 
various charge and/or spin distributions in low temperatures. 



II. METHOD OF CALCULATION 



In this paper we consider only the case of p/ = 1/2 and = 1. Then, for the system composed of N sites there are 
[(A72)!]^ possible charge distributions and, for each of them, there exist 2^^^ spin configurations. So the entropy (per 
site) resulting from various possible distributions of f-electrons is a sum of the charge part equal to j^Log ^(■jy^J)! ] ^ ' 
and the magnetic part equal to ^Log2. In the limit of — > 00 the charge part of entropy tends to Log2. 

For any fixed configuration of localized electrons C one can calculate all thermodynamic characteristics of the 
d-electron subsystem from standard formulas (see for example Q). In particular, the Gibbs potential per site gc is 
given by the expression 

5^(^) = -]^E^«5(l + e-''[^"(^)--l) (2) 

n 

and the free energy fc by 

fc{Pd) = gciKPd)) + lJ'{Pd)pd (3) 
with the following condition fulfilled (in our case pd — 1): 

ly \ (4) 



Pd = 



Then, the partition function per site z is obtained from the summation over all configurations C of f-electrons 



and the free energy / is equal to 



^E^-^^^" (5) 
c 



f--^LogiY^e-^fo-). (6) 



Since we are not able to determine / exactly in the whole temperature range, we replace it by three different approx- 
imate functions, relevant for three different phases: CSO, CO and NO. In order to do this we first divide the sum in 
JS]) as follows 

^e-^^/^ =e-^^^^sco + J2 ^'^"^^ + E e^'"'^^^. (7) 
c Cdc:^ cec- 

where fcsco ^^e free energy of the SCO phase, C*^ denotes the set of checkerboard configurations Cch (with respect 
to charge) of the f-electrons with all possible spin arrangements but SCO and C* denotes all C configurations but 

Cch- 
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FIG. 2: Free energies fc (calculated per site) of phases with fixed configurations of the f-electrons (the continous lines) 
corresponding to the highes and the lowest value functions taken from a choosen set of configurations: a) C G Cch b) all C. 
The points mark approximative average free energy function taken in the calculations: a) (/c)csCchi b) {fc)auc- 



With an increase of temperature the subsequent terms in ([7]) contribute to /: for the lowest temperatures only the 
first term, then, for intermediate temperatures two first terms and, finally, for higher temperatures all three terms. 
So for temperatures very close to zero the main contribution to / comes from the free energy fccso^ attributed to 
the ground state configuration CSO of the f-electrons 

fc^-^^Log{e-f'f^cso)^fc^^^. (8) 

For higher temperatures we calculate mean free energies by averaging them over appriopriate sets of J- electron con- 
figurations. In order to find the mean free energies we take into account a representative class of configurations C 
of the f-electrons only, composed of all periodic phases with periods not exceeded four lattice sites. Then our trial 
set contains 13 non-equivalent configurations. Three of them belong to the Cch class. For all these configurations 
we calculated exact free energies per site in the limit of large N by solving the eigenvalue problem and finding the 
eigenvalues £"„ (C) . This required us to diagonalize up to 4 x 4 matrices and result in analytical formulae for at most 
4 different energy bands, separately for spin-up and spin-down d-electrons (for more details see Refs. 2, 5, 6]). 

For intermediate temperatures we tak an average fc over the family Cch only. Their free energies fc are differrent, 
but the difference is quite small, even for T = 0, and tends to zero with temperature (see Fig. 2). 

The whole contribution to / coming from this family could be expressed by a mean value {fc)cGCah the following 
way. 

/ = /co = -^Lo5(2^/^e-''^</->---c. ) = -^Log2 + {fc)cec., (9) 

We get the appropriate mean free energy {fc)ceCah using the arithmetical average of the relevant free energies fc- The 
averaging procedure we used is equivalent to taking the geometrical mean value of relevant ingradients of the partition 
function. Indeed, we get (O if we replace each element of the sum ^ g-P^fc j^y Mc^^ ■ ( H e^^^^'^y^^'^'^'<!h ^ 

where Mc^^ = 2^/^ is the number of configurations in Cch 



4 



The same procedure we use when calculating an approxiamtive mean free energy /ato of the non-ordered phase. 
However, then we take an arithmetric average of free energies of all configurations from the choosen set, and take into 
account the fact, that in the large N limit there are 2^/^ [(^)!]'' ~ 2^/^2^ such configurations. So we have 

f = fNO = -^Log2 - ^Log2 + {fc)auc. (10) 

We tested the averaging procedure for the spinless Falicov-Kimball model and found, that the transition temperatures 
we got are nearly the same as those obtained by the Monte Carlo simulations in the whole range of the interaction 
parameter U Q. 

III. RESULTS AND CONCLUSIONS 

The results of our calculations are displayed in Fig. 3. In the upper panel a) it is shown the magnetic order-disorder 
transition temperature Tmo between SCO and CO phases and in the lower panel b) the charge order-disorder transition 
temperature Tco between the CO and NO phases. Both Tmo and Tco were calculated for J = 0.2U as functions of 
the rescaled on-site interation parameter jj^- 

. 004 
. 0035 

. 003 
. 0025 
i 0.002 
. 0015 

0. 001 
. 0005 

0.6 0.7 0.8 0.9 
U/ (U + 1) 

0.3 
0.25 

o 

u 0.2 

H 

0.15 
. 1 

0.6 0.7 0.8 0.9 1 
U/ (U + 1) 

FIG. 3: The temperatures of order-disorder phase transitions; a) magnetic Tmo and b) charge Tco, as a function of the on-site 
interaction parameter (in the both cases J — 0.2(7). 

Since we used rather crude approximative scheme of calculations of the free energies, the obtained values of the 
transition temperatures could not be very precise. However, we expect that the quantitative picture is correct. In 
particular, two main conclusions seem to be relevant. First, the maxima of the both transition temperatures are 
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attained at very close to each other values of the parameter w 0.8 -j- 0.85. It corresponds to J/ » 4 5.67 (and 
J~0.8-^1.13, respectively). 

Another important message resultant from the calculations is the fact, that Tmo appears to be almost two orders 
of magnitude lower than Tco, even though the Hund coupling constant J is assumed to be only 5 times smaller 
than U. The possible cause of this rather unexpected outcome is the lack of direct interaction between the spin-up 
and spin-down itinerant electrons (the Hubbard term) in the model. We anticipate that the inclusion of this type of 
interaction would enhance the magnetic couplings in the system, resulting in an increase of the magnetic transition 
temperature. 
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